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The  scattering  of  plane  congressional  waves  by  a  spherical  obstacle  in  an 
elastic  solid,  which  was  investigated  by  Yii«  and  Truell  is  further.  For 

a  rigid  inclusion,  the  boundary  conditions  are  redefined  to  take  into 
the  motion  of  die  inclusion  inside  die  solid.  By  a  proper  Hmi+iwy  process,  it  is 
shown  that  the  solutions  for  a  rigid  insert,  a  fluid  sphere,  a  cavity,  or  an  obstacle 
in  a  fluid,  are  all  derivable  from  the  general  results  of  an  elastic  iw«inT«"n,  in 
each  case,  the  rates  of  energy  scattering  due  to  a  small  obstacle  are  found  to  be 
inversely  proportional  to  the  fourth  power  of  wave  length. 
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ABSTRACT 

The  scattering  of  plane  compressional  waves  by  a  spherical  obstacle  in  an 
elastic  solid,  which  was  investigated  by  Ying  and  Truell  is  examined  further.  For 
a  rigid  inclusion,  the  boundary  conditions  are  redefined  to  take  into  consideration 
the  motion  of  the  inclusion  inside  the  solid.  By  a  proper  limiting  process,  it  is 
shown  that  the  solutions  for  a  rigid  insert,  a  fluid  sphere,  a  cavity,  or  an  obstacle 
in  a  fluid,  are  all  derivable  from  the  general  results  of  an  elastic  inclusion.  In 
each  case,  the  rates  of  energy  scattering  due  to  a  small  obstacle  are  found  to  be 
inversely  proportional  to  the  fourth  power  of  wave  length. 
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I.  INTRODUCTION 

The  scattering  of  a  plane  compressional  wave  by  a  spherical  obstacle  in 
an  elastic  solid  infinite  in  extent  has  been  investigated  by  Ying  and  Truell.  “ 

Three  types  of  obstacles,  an  elastic  sphere,  a  rigid  sphere,  and  a  cavity  are 
discussed;  detailed  results  being  given  for  the  case  of  Rayleigh  scattering,  i.  e. , 
when  the  size  of  the  obstacle  is  much  smaller  than  the  incident  wave  length.  The 
solution  of  the  elastic  sphere  is  quite  general  in  the  sense  that  it  applies  for 
various  values  of  material  constants  (Lame' s  constants,  density  etc.)  character¬ 
izing  the  scatterer  and  the  surrounding  medium.  However,  Ying  &  Truell 
studied  the  problems  of  the  rigid  sphere  and  a  cavity  separately  because  some  terms 
in  the  general  solution  are  ambiguous  when  the  limiting  values  of  the  density 

of  the  inclusion  are  taken.  For  the  same  reason,  the  fourth  type  of  obstacle,  a 

2 

cavity  filled  with  an  inviscid  fluid,  was  also  treated  independently. 

Moreover,  the  scattering  by  a  rigid  sphere  was  calculated  with  the  assumed 
boundary  condition  that  the  displacement  due  to  the  combined  Incident  and  re¬ 
flected  waves  vanishes  at  the  surface  of  the  sphere.  This  led  to  the  conclusion 
that  the  energy  scattered  by  a  small  rigid  sphere  per  unit  time  is  independent 
of  the  incident  wave  length;  this  result  being  quite  different  from  other  well 
known  scattering  phenomena,  i.  e. ,  the  rate  of  energy  scattering  is  inversely 
proportional  to  the  fourth  power  of  wave  length. 

The  same  boundary  conditions  for  a  rigid  obstacle  have  been  applied  in 
3  4  5 

many  related  investigations.  ’  ’  In  a  recent  study  on  the  dynamical  stress 

0 

concentration  in  an  elastic  plate  with  a  circular  rigid  insert,  the  authors  find 
that  if  the  zero  displacement  boundary  conditions  are  used,  stresses  will  become 
infinitely  large  at  certain  points  in  the  plate,  a  conclusion  even  more  drastic 
than  the  unusual  rate  of  energy  propagation  mentioned  above. 
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In  a  fluid  medium,  the  disturbance  of  plane  waves  (sound  waves)  due  to  a 
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rigid  sphere  was  first  treated  by  Rayleigh.  He  showed  that  for  a  fixed  sphere, 
there  was  a  net  unbalanced  force  acting  on  the  surface  of  the  sphere,  i.  e. ,  if 
the  sphere  were  not  held  fixed  by  some  external  constraint,  it  would  move.  A 

g 

modified  solution  for  an  unconstrained  rigid  sphere  was  given  by  Lamb. 

However,  in  either  case,  the  rate  of  energy  scattering  is  inversely  proportional 
to  the  fourth  power  of  wavelength. 

There  is,  however,  a  major  difference  between  a  movable  and  a  fixed 
obstacle  in  an  elastic  solid.  With  the  view  of  obtaining  a  mechanical  Illustration 
of  the  selective  absorption  of  light  by  a  gas.  Lamb  also  investigated  the 
scattering  of  shear  waves  in  an  incompressible  solid.  If  the  boundary  of  the 
sphere  is  absolutely  fixed,  the  rate  of  energy  scattering,  called  by  Lamb  the 
"dissipation  ratio, "  is  independent  of  wavelength.  Lamb  explained  that  this  was 
accounted  for  by  the  abnormal  degree  of  constraint  imposed  on  both  radial  and 
tangential  displacements  of  the  surrounding  medium,  whereas  in  the  corresponding 
fluid  problem,  there  was  complete  freedom  of  lateral  motion  at  the  surface  of 
the  rigid  sphere.  When  the  sphere  is  allowed  to  move  with  the  surrounding  medium, 

9 

the  dissipation  ratio  conforms  with  the  general  rule.  Later,  Sezawa  investigated 
the  problem  of  scattering  of  elastic  waves. 

Nishlmura  and  Jimbo11  also  investigated  the  diffraction  of  plane  com- 
presslonal  waves  by  three  types  of  scatterer,  l.e. ,  elastic,  rigid  and  cavity. 
However,  as  a  result  of  the  assumed  standing  waves,  only  even  terms  of  series 
expression  for  waves  in  polar  coordinates  are  needed,  thus  no  complication  such 
as  that  encountered  in  the  case  of  traveling  waves  arises. 

In  this  paper,  the  boundary  conditions  for  rigid  scatterer  are  redefined 

12 

in  order  to  take  into  account  the  rigid  body  translation.  It  is  then  shown  that 
by  a  limiting  process,  the  results  for  a  rigid  sphere,  a  cavity  or  even  a  fluid 
sphere  can  all  be  derived  from  the  general  solution  of  scattering  due  to  an  elastic 
inclusion.  The  process  can  be  applied  further  to  deduce  the  results  of  the 
scattering  of  sound  in  fluid.  The  rates  of  energy  scattering  are  found,  in  all 


cases,  to  be  inversely  proportional  to  the  fourth  power  of  the  incident  wave 
length. 
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In  the  interests  of  brevity,  only  results  which  are  to  be  compared  with  that 
in  References  1 ,  2  or  which  are  essential  to  carry  out  the  limiting  process  are 
given.  For  a  detailed  discussion  and  references,  the  reader  is  referred  to  the 
original  work  by  Ying  and  Truell. 

n.  BASIC  EQUATIONS 

In  an  elastic  solid  with  waves  propagating  symmetrically  about  the  z  axis 
(Figure  1),  the  displacement  vector  in  spherical  coordinates  (r,  0,  <P)  can  be 
expressed  in  terms  of  two  displacement  potentials  4  and  * ,  i.  e. , 

u  =  V*  +  Vx  (e^  d’t/dd)  (1) 

where  V  is  the  gradient  operator  and  is  a  unit  vector  along  <P  coordinate 
curve.  Each  potential  can  then  be  shown  to  satisfy  a  scalar  wave  equation 

c^V2*  =  82* /at2  (2) 

c2  V2  *  =  92*/3t2  (3) 

r  1/2  1/2 

where  =  [  (^  +  2p)/ p  J  ,  c^=  (p/p)  are  the  velocities  of  compressions: 

(longitudinal)  waves  and  shear  (transverse)  waves  respectively;  A,  p  are  Lame '  s 
constants,  and  p  the  density  of  the  solid. 

The  stress  components  are  related  to  the  displacements,  in  cfyadic  notation, 
as 

r  =  A(V  •  u)i  +  2p(Vu  +uV)  (4) 

where  t  ,  1  represent  the  stress  tensor  and  isotropic  tensor  respectively.  Due 
to  axial  symmetry,  the  displacement  u^  as  well  as  the  stresses  T T(^  >  Tg^ 
vanish. 
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When  plane  waves  impinge  on  the  surface  of  an  inclusion  embedded  in  an 
infinitely  extended  elastic  solid,  two  types  of  waves  (compressional  and  shear) 
are  reflected  back  into  the  infinite  medium  and  two  waves  are  refracted  into 
the  inclusion.  For  the  convenience  of  ensuing  discussion,  the  infinite  solid  is 
designated  as  medium  1  and  the  spherical  inclusion  as  medium  2.  All  material 
constants  such  as  X,  M,  P  etc.  in  each  medium  will  be  distinguished  by  subscript 
1  or  2  accordingly.  The  potentials,  displacements  and  stresses  associated  with 
the  incident  waves  will  be  designated  by  the  superscript  (i);  those  with  the  re¬ 
flected  waves  by  (r)  and  those  with  the  refracted  waves  by  (f) . 


m.  INCIDENT,  REFLECTED  AND  REFRACTED  WAVES 


Let  the  incident  plane  compressional  waves  traveling  in  the  positive  z 
direction  be  represented  by  two  potentials  (Figure  1 ) 


j  .«  -  vi(V-“*> 

f  =  0 

where  (=  2ir/  wave  length)  is  the  wave  number  of  compressional  waves  in 
medium  1,  and  w  is  circular  frequency.  In  order  to  be  a  solution  of  the  wave 
Equation  (2),  the  wave  number  has  the  following  value 


“l  =  <j2pl/  +2/1i>  • 

2 

From  (1)  it  is  clear  that  *  ,  a  constant,  has  the  dimension  of  (length)  and 

^  13 

a  *  is  the  amplitude  of  the  incident  wave.  In  spherical  coordinates, 


(2n  +  l)i“j  (a  r)P  (cos  0), 
n  l  n 


(6) 
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In  which 


V*  '  ffx  J„  +  1/2(X> 

is  the  spherical  Bessel  function  of  the  first  kind,  P  (x)  is  the  Legendre  polynomials. 

-i«t  n 

In  Equation  (6)  and  sequel  the  time  factor  e  is  omitted  whenever  its  occurrence 
is  apparent.  By  substituting  (5)  or  (6)  into  (1)  and  (4),  the  displacements  and 
stresses  of  the  incident  waves  can  be  determined. 

Expressions  for  the  reflected  and  refracted  waves  can  be  obtained  from  the 
solutions  of  the  wave  Equations  (2),  (3).  In  an  infinite  solid,  only  waves  propa¬ 
gating  outward  from  the  center  will  be  considered.  Thus  the  reflected  waves 
are  given  by 

•0 

*(r)  =  ^  Anhn(air)  Pn(cos  0)  (7) 

n=0 

•0 

*<r)  -  X  Bn  hn  <V>  Pr> <C°"  S)  (fl 

n=0 

with 

hn(x)  =  hn(1)  (X)  -  Hn+1/2(1,(X) 

being  the  spherical  Hankel  function  of  first  kind.  As  there  is  no  confusion  with 
the  Hankel  function  of  second  kind  which  represents  converging  wave  in  this 
study,  the  conventional  superscript  (1)  is  omitted. 

The  refracted  waves,  being  confined  in  the  spherical  scatterer,  are 
standing  waves.  They  can  be  represented  by 


(9) 


6 

00 

*(1)  -  - 1  cni»c  v>  p„<COB  •> 

n=0 


*<f>  -  -  I  DnJ..C^r>  P„  (oos  *> 

n=0 


(10) 


In  (7)  -(10), 


<*\  =  w2pl/(Xi  +  S>ii) 

1-1,2  (11) 

02  2  / 

and  A  ,  B  ,  C  ,  D  are  coefficients  to  be  determined  from  the  boundary  condi  - 
n  n  n  n  J 

tions  at  the  surface  of  the  inclusion. 


In  the  medium  1,  the  resultant  waves  are  then  determined  by  superposing 
the  incident  and  reflected  waves  whereas  fclie  refracted  waves  are  the  only  ones 
in  medium  2.  Displacements  and  stresses  in  each  medium  are  summarized 
below. 


*  =  *(*)  *(r>  • 
f  ^  ^  (r) 


(12) 


u  =  -  )  r-*A«,+A  K,  +  B  P  (cos  6) 

rl  r  L/  L  0  1  n  11  n  12  J  n'  ' 

n=0 

(13) 

1  V  r  *i  dP_(cos  9) 

UM=?Z  [-*0  52  +  An5321+Bn  522]  - di -  * 

n=0 
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I 

I 

I 

I 

I 

l 


7- 


1 


I 

I 

I 

1 

I 

I 

I 

I 

I 


rrl  2 
r 


2Mj  ^ 

~J  L  [‘*0  VAo*Sl  + V32]  p„<00*,>  ’ 


n=0 


r01 


2/ii  V  r  dP  (ooi 

=  2  /  r-*n5.+A  S.i+B 

2  O  L  0  4  n  41  n  42  J  d0 

n=0 


(cos  0) 


(14) 


$  = 


*  _ 
2 


(15) 


Ur2  =  '7  I  [CnS13  +  Dn*14]P„<M“',> 


n=0 


i  r'  „  of  i 

u  =  -  —  /  r  c  &  +  d  &  i  — - 

32  r  L  n  23  n  24  J 


dP  (cob  0) 


d0 


n=0 


rr2 


2^2  V 

”TZ  LCn«33  +  D„‘34]  Pn<C°*  »> 


n=0 


(16) 


(17) 


,  .  -2i 

r  02  2 

r 


Y'  ay  i 

/  r°  + °  s,,.  l  — - 

Li  L  n  43  n  44  J 


dP  (cob  0) 


dO 


n=0 


where 
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*1  ■  -i“<2,1+1)[»i„(“1r)  "“lrin  +  l(“lr)] 

*2  “  "  l”  (2n  +  1)  )„(“!*) 

*3  ‘  _  l“  (2“  +  '>  [<“2 -f  yoj*)  tSOjr  )n+1(“jr)] 

*4  -  -  l"(2n  +  l)[(n  -  1)  yiy)  -  'V'n  +  ^V*]  ' 

«n  =  »v  V>  '  V  WV> 

*21  =  W> 

*31  *  <”2-n-|^il2)  VV,+!VW01I> 

*41  =  <n  -  1)  VV>  -  VVi<V>  • 

S12  -  -  n(n  +  1)  hn(^1r) 

*22  =  -(n  +  1>h„(V,+ 

*42  ■  -C‘-1-|^/>VV»'V,|l-U<V)  ' 

*13  =  “W  '  V  WV> 

*23  * 

*33  =  <"2  '  \  IA1  +  2VWV> 

*43  ■ (n  -  *>  vv>  -  V  WV>  • 

S14  =  *■  n(n  +  1)  jn(/?2r) 

*24  = 

*34  ■ 

*44  =  -  C2  -  1  -  i  ln^r)  -  VWV>  ■ 


9 


Note  that  5  ^  =  <Re  ( 5  tl)  and  5  ^  =  <Re  (  £  provided  the  subscripts  of  wave 

numbers  o^,  are  changed  prbperly.  All  £  1  s  are  dimensionless  numbers 

and  A  etc.  have  the  same  dimensions  as 

n  0 


IV.  ELASTIC  INCLUSION 


If  the  spherical  inclusion  of  radius  a,  is  elastic,  and  it  is  bound  to  the 
surrounding  medium,  the  tractions  and  displacements  must  be  continuous  at 
the  interface.  Thus  at  r  =  a,  the  continuity  conditions  are 


(i)  ^  (r)  (f) 

u  +  u  =  u 


(i)  A  (r)  (f) 

uV + uV  ■ u  V 

t<‘>  +  r<r>  .  r<f> 

rr  rr  rr 

T(i)  .  T(r)  _  T(f) 

re  re  re 


(19) 


Since  the  Legendre  polynomials  Pn(cos  e)  and  the  associated  Legendre  polynomials 


Pn  (C°8  6)  =  '  dO  Pn(C0S  0) 


each  form  an  orthogonal  complete  set,  when  (13),  (14),  (16)  and  (17)  are  sub¬ 
stituted  into  Equations  (19),  the  coefficients  of  the  Legendre  polynomials  on  both 
sides  of  (19)  must  be  equal  for  each  value  of  n.  This  results  in  four  simultaneous 

algebraic  equations,  sufficing  to  determine  the  four  unknown  coefficients  A  ,  B  , 

•  n  n 

C  ,  and  D  .  In  matrix  form,  these  equations  are 
n  n 


E 

E 

E 

E 


11 

E12 

E13 

21 

E22 

E23 

31 

E32 

pE33 

41 

E42 

PE43 

E14 

A 

n 

E24 

B 

n 

PE34 

C 

n 

PE44- 

D 

L  n  J 

(20) 
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where 


and 


E«  =  (S«’r  =  a 


P  -  V"! 


(21) 


Except  for  the  limiting  cases  when  the  arguments  of  the  spherical  Bessel  functions 
become  very  small  or  very  large,  these  equations  can  best  be  solved  numerically. 

The  scattering  of  waves  by  small  obstacles  (a^a  «  1)  is  generally  known 
as  Rayleigh  scattering.  From  (11)  it  follows  that  if  a  a  « 1,  a  &  and  /3  a 

±  a  1 

will  all  be  very  small.  Thus  the  spherical  Bessel  and  Hankel  functions  in  (20) 
can  be  approximated  by  the  leading  terms  of  their  series  expansions  which  are 


1  <x)  =  (2x)° 


V  (-l)r  (n  +  rl  !  2r 

Z  r !  (2n  +  2r  +  1)  !  x 

r=0 


ix 

h  (x)  =  -g—  > 

n'  '  n  +  1  Lt 
x 

r=0 


(-i)n~r  +  1  (n  +  r)  !  n-r 
r!  (n-r)  !  2r  X 
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Thus  for  the  case  of  Rayleigh  scattering, 


A0  = 


Ji* 

3  0 


1  - 


L  8p(^2/«2)  -4(P-1)J 


(“ja)3, 


A1  =  3  *0  (1  "  P2/Pl)  <ai“)3' 

VV> 


2n-l 


n  >1, 


(22) 


N  = 


(4n  -l)(p-l) 


n[<"+1>p+ 


P  -  1  + 


B  =  -  (0  /a  )n  +  1  A  /n 
n  11  n 


n  >0 


The  Coefficient  is  not  needed  in  the  calculation  of  stress  and  displacements. 

If  is  set  equal  to  (-1)  i/o^,  the  above  results  agree  with  those  of  Reference  1. 
It  is  worth  to  note  here  that 


V2*i 


2(1  -  a) 
1  -  2cr , 


i  =  l,  2 


where  c  are  the  Poissons'  s  ratios  of  the  media. 
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V.  FLUID  INCLUSION 

An  inviscid  fluid  will  neither  resist  shearing  stress  nor  support  shear 
waves.  If  a  cavity  in  an  elastic  solid  is  filled  with  such  a  fluid,  only  compressional 
waves  are  refracted  into  the  fluid.  The  boundary  conditions  at  r  =  a  are  reduced 
from  (19)  to 


U<1>  + 

(r)  (f) 

u  =  u 

r 

r  r 

r<i>  + 

T(r)  =  T(f) 

rr 

rr  rr 

T<i>  + 

T&  =  o 

rd 

r6 

The  circumferential  displacement  u  in  the  solid  might  be  different  from  that 
in  the  fluid  at  the  interface.  Such  a  paradox  in  discontinuity  is  due  to  the  assumed 
zero  viscosity. 

The  displacements  and  stresses  in  the  elastic  solid  are  obviously  still  the 
same  as  (13),  (14)  with  Aq,  satisfying  (23).  In  the  fluid,  they  can  be  obtained 
from  (16),  (17)  by  taking  the  limiting  values  as  m2  °- 

It  should  be  noted  from  (11)  that  as  M2  “*  0,  «  ;  but 


remains  finite. 


with 


„2  2  a2  , 


(24) 


By  setting  ju  =  0,  D  =0  (no  refracted  shear  wave),  and  taking 
m  n 


M2  S33  M1  S33  ’ 


5  33  =  ’  J 


(25) 
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one  finds  that  the  boundary  conditions  (23)  require  that  Aq,  Bq,  satisfy 
three  equations,  i.  e. , 


E11  E12  E13 


E  E  E 
31  32  33 


EE  0 
41  42 


A 

X ' 

n 

1 

B 

= 

E„ 

n 

0 

3 

C 

J 

.  n  . 

L  4j 

(26) 


2 

These  equations  agree  with  that  given  by  Einspruch  and  Truell. 

The  same  equations  can  directly  be  deduced  from  the  case  of  an  elastic 
inclusion.  Recall  that  B—°°  as  u  —  0;  hence  the  spherical  Bessel  func- 

A  « 

tions  with  argument  /3  a  in  the  matrix  (20)  can  be  approximated  by  the  asymptotic 
formula 


j  (x) 

n  ' 


X"*  * 


1 

X 


cos  (x  - 


n  +  1 
2 


it)  . 


It  is  clear  that  as  —  • , 


yv> =  o(e> 

Vjn  +  l(^2a)  =  °(e<)) 

*2  =  o(e 2) 


where  e  is  an  infinitesimal  quantity  and  ”o  (e)  "  means  "the  order  of  magnitude 
of  e  Applying  these  limits  to  all  elements  of  the  matrix,  one  finds  that 


=  0(€), 

pE44  =  o(e)  , 


pE 


34 


pE 


43 


pE 


33 


o(e2)  , 
o(e2)  , 
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By  neglecting  all  elements  with  order  of  magnitude  higher  than  ,  Equations 
(20)  reduce  to 


E,  ~  E„ .  E.  *  0 

'  A 

11  12  13 

n 

0 

E  E  E  o(c) 

21  22  33  '  1 

B 

n 

E_,  E„„  E*  0 

C 

31  32  33 

n 

E  E  .  0  0 

D 

9 

i 

L  a  J 

(27) 


For  the  determination  of  the  coefficients  A  ,  B  ,  and  C  ,  these  equations  can 

n  n  n 

easily  be  reduced  to  (26). 

For  Rayleigh  scattering,  A^,  can  be  directly  deduced  from  (22)  by 
taking  the  limiting  values  asp— 0,  P  —  ».  They  have  the  following  values 

a 


which  check  with  those  given  by  Einspruch  and  Truell. 

2 


2,14 


A0  -  - 1  i*0 


3P2P\/Pia2  +  4 


A1  =  3  *0  <1‘p2/pi)(aia)3’ 


n  >1 


V 


4n2  -  1 


2  /  B  \2  » 

_  2n  +1  ri_\ 

n(2n  -2)  y  ) 


(28) 


n  >0  . 
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VI.  RIGID  SPHERE 
1.  Boundary  Conditions 

When  the  inclusion  is  much  more  rigid  than  the  surrounding  material,  it 
may  be  treated  as  a  perfectly  rigid  sphere.  By  definition,  the  distance  between 
any  two  points  in  a  rigid  body  remains  constant  at  all  times,  but  because  of  the 
surrounding  medium  being  elastic,  the  sphere  will  translate  as  a  rigid  body  under 
the  impact  of  the  incident  waves.  Unless  the  sphere  is  fixed  in  position  by 
external  forces  or  other  means,  which  is  hardly  conceivable  when  it  is  embedded 
in  an  infinite  solid,  the  displacements  of  the  elastic  medium  should  be  equal  to 
the  translation  of  the  inclusion  at  the  interface. 

Let  U  denote  the  rigid  body  motion  of  the  sphere  along  the  direction  of 
z 

incident  waves  (z  axis).  The  boundary  conditions  at  r  =  a  become 


U  cos  0  =  U  P  (cos  0) 

Z  Z  1 


.<*>  +  U(r> 


e 


e 


-  U  sin  0  = 
z 


U 


dP^cos  0) 


d0 


(29) 


The  translation  U  is  to  be  determined  from  the  equation  of  motion 

Z 

mUz  =  ^  (t^cob  0  -  Tr0sin  0)*^  sin  0  dfld  <t>  (30) 

4  3 

where  m  =  r  ira  p  being  the  mass  of  the  sphere  and  the  integral  over  the 
spherical  surface  representing  the  force  component  acted  on  the  sphere  by  the 
surrounding  medium. 
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15 

Restoring  the  time  factor  and  substituting  (14)  into  (30) ,  one  finds 

Vf  [31  V!<V>  +  Alhl(“la)  "  2B!  W>]  e  iWt  (31) 

with 

rj  =  P/P2  • 

Owing  to  the  orthogonality  of  P  (cos  6)  and  P 1  (cos  0),  the  terms  with  n  #  1  in 

n  n 

the  series  (14)  vanish  after  integration. 

The  boundary  conditions  used  in  Reference  1  and  others  are  special  cases 
of  (29)  with  rj  =  0,  which  implies  that  the  density  of  the  sphere  is  infinitely  large. 
However  no  consideration  of  rigid  body  motion  is  needed  if  only  standing  waves 
are  involved  as  in  Nishimura  and  Jimbo1  s  work. 11  In  that  case,  the  sphere, 
even  with  finite  density,  is  always  being  held  in  the  equilibrium  position  by  two 
trains  of  waves  propagating  in  opposite  directions. 

2.  Solution 

For  n  #  1,  the  coefficients  A  ,  B  are  determined  from  the  matrix 

n  n 

equations 


Coefficients  are  given  by  a  different  set  of  equations 
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with 

en  -  (1  -  rfl  h^o^a)  -  a1ah2(a1a) 
e21  =  (1  -  tj)  h1(a1a) 
e12  =  -2(1  -rjjh^a) 

(34) 

®22  =  -2(1 -TjJh^aJ  +  ^ah^a) 
ex  =  -  31  [  (1  -  77)  j^a)  -  a1aj2(“1a)] 
e2  =  -  3i  (1  -  tj)  J1(«1a) 

Ifthescatterer  is  very  small  in  comparison  with  the  incident  wave  length, 
the  coefficients  are  found  as 


-J‘VV>3 


A,  = 


A  =  i 


I  *0  (1  -  P2/Pt)<“1a»3 

r  n  "1 2 

n-1  .  2  n!  ..  .  .  2n-l 

#o[WJ  W> 


N  = 


n(4n_T  1) 


n  +  (n  +  1)  (^/aj ' 


n  >1 


(35) 


B  =  -  (P. /a  ) n+1  A  /n  n  >0 

n  'll  n 


A  and  B  are  identical  with  those  given  in  Reference  1  except  for  n  =  1. 
n  n 
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3.  Derivation  from  Solutions  of  Elastic  Sphere 

The  lengthy  calculations  of  the  rigid  body  translation  given  in  (31)  can  be  * 

avoided  if  the  solution  of  an  elastic  inclusion  is  known  in  advance.  Caution 

should  be  exercised  however  in  deriving  the  results  from  Section  IV.  Although 

it  is  well  known  that  no  waves  propagate  in  a  rigid  solid,  one  cannot  simply  set 

all  C  ,  D  of  the  refracted  waves  to  be  zero  in  this  study, 
n  n 

In  (16)  or  (13)  for  constant  r,  there  are  two  terms  in  the  series  having  the 

form 


uf  =  bjP^cos  B)  =  b^  cos  Q 
ue  .  -cjdp 


(cos  0)/dfll  =  c  sin  B 


'  J 


(36) 


where  b^,  are  constants.  By  transformation,  the  corresponding  displacements 
in  z  -  q  coordinates  (Figure  1 )  can  be  shown  as 


u 

z 


+  c1)  cos  2  6 


Uq  =  !<bl  +  Cl>  Sin2fl 

There  is  clearly  a  rigid  translation  ~  (b  -  c  )  in  the  z  -  direction  implied  in 

M  1  i. 

(36).  Indeed,  this  rigid  body  translation  can  be  separated  from  the  deformation 

16 

if  (36)  is  rewritten  in  the  following  form 


u 

r 


\  (bx  +  Cj)  cos  6  +  |  (bx  -  cx)  cos  B 


u 


\  (bx  +  Cj)  sin  8  -  |  (bx  -  c^  sin  6 


(36a) 
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By  suppressing  all  terms  (including  n  =  1 )  in  the  series  expressions  for  u^p  , 
u<f> ,  the  sphere  is  then  assumed  to  be  rigid  and  fixed  in  space.  For  a  movable 
sphere,  this  simple  suppression  will  lead  to  erroneous  result.  In  fact,  the 
second  terms  in  (36a)  are  exactly  what  have  been  incorporated  in  setting  up  the 
boundary  conditions  (29). 

A  limiting  process  analogous  to  that  in  Section  V  will  now  be  used  to 
derive  the  rigid  sphere  solution.  If  the  medium  2  is  perfectly  rigid,  it  is  seen 
that 


but  p  remains  finite.  Thus  by  taking  only  the  leading  terms  of  the  series  of 
spherical  Bessel  functions,  one  finds 

a2a  =  o(e),  /^a  =  o(e), 

p  =  o(e  ), 

W  =  °(eI1)’  V3n+l(^2a)  =  0(cI1+2>- 

If  these  limiting  values  are  substituted  into  the  matrix  (20) ,  it  can  be 
shown  that  for  n  #  1; 

*U  -  °«n>.  1-1.  2;  3,  4 

pSy  =0(£n-2),  k,  «  =  3,4 

for  n  =  1, 

Eil  =  0(€)’ 
pEkl  =  °(c)’ 

It  follows  then,  for  n  #  1,  the  matrix  reduces  to  (32)  when  are  neglected 
in  comparison  with  pE^  .  For  n  =  1 ,  the  complete  4x4  matrix  should  be  used 
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with  the  spherical  Bessel  functions  in  the  third  and  fourth  columns  being  replaced 
by  the  corresponding  small  argument  approximations.  Coefficients  A1 ,  thus 
determined  are  found  to  be  the  same  as  that  from  Equations  (33).  Furthermore, 
all  coefficients  for  the  case  of  Rayleigh  scattering  as  given  in  (35)  can  directly 
be  deduced  from  (22)  by  letting  p  approach  infinity. 

Although  in  this  limiting  process,  C  ,  D  are  treated  as  finite  numbers, 

n  n 

the  refracted  waves  still  vanish  because  of  the  small  values  for  /3  a,  a  a  which 

2  2 

appear  in  the  series  representations  of  $  ,  't  . 

2  2 

Vn.  SPHERICAL  CAVITY 

Unlike  a  rigid  inclusion,  the  case  of  spherical  cavity  presents  no  compli¬ 
cation  in  analysis.  Coefficients  A^,  B^are  fixed  by  the  boundary  conditions  at 
r  =  a, 

T(l)  ♦  .  0 

rr  rr 

(37) 

♦  P-t, 

rO 

which  yields  two  simultaneous  equations 


E 

31 

E32 

— 

> 

p 

i — 

n 

W 

1 _ 

E 

41 

E42_ 

B 

n 

o 

e 

ii 

- 1 

x* 

EC 

_ 1 

(38) 


These  equations  can  also  be  derived  from  (20)  by  a  limiting  process  —  0, 

H  —  0,  p  —  0)  similar  to  that  used  in  Section  V. 

2  2 

With  small  a^a,  Aq  and  Bq  can  be  obtained  by  either  solving  Equation  (38) 

or  by  taking  p  — ■  0,  p  -*  0  in  (22).  They  are  listed  below  for  reference. 

2 
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A 

n 


,n 


N  (a  a) 

p  1 


2n-l 


n  >1 


N 


P 


4n  -  1 


1  - 


2n2  +  1 
2n(n  -  1 ) 


B 

n 


-  (0/^) 


n+1 


A  /n 
n 


n  >0 


(38) 


Vffl.  SCATTERER  IN  FLUID 


By  applying  an  analogous  limiting  process  of  reducing  the  shear  rigidity 
of  medium  1  (Mj  —  0)  to  the  solutions  given  hitherto,  the  scattering  of  sound 
(compressional)  waves  in  a  nonviscous,  compressible  fluid  can  be  obtained. 

The  obstacle  could  be  an  elastic  sphere,  a  rigid  sphere  (fixed  or  movable),  a 
cavity  (bubble  in  liquid)  or  an  acqueous  sphere  (hydrosol).  Most  of  these 
solutions  have  long  been  in  existence.  It  is  sufficient  here,  as  an  example,  only 

to  discuss  the  scattering  of  sound  by  a  rigid  sphere,  which  was  first  treated 

7  8 

by  Rayleigh  and  Lamb. 


Because  in  the  inviscid  fluid,  no  shear  wave  can  propagate,  the  shear 
modulus  — '0  and  the  shear  wave  number  ®  .  By  the  application  of 

the  asymptotic  formula  of  spherical  Hankel  functions  for  large  argument 


h  (x) 
n 


(-i) 


n  +  1  e 


ix 


00 


X 
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It  can  be  shown  that,  in  (32) , 

V  =  o(€_1),  hn(^a)  =  o(e); 

E12  =  °<€> 

E22  '  °<e<>> 

If  E  is  neglected  in  comparison  with  Eno,  Equation  (32)  degenerates  into, 


A.  E 
n  11 


-  Vi 


(40) 


and 


A 


n 


-♦0l”(2n-H) 


ny  v1  -  vwv> 

nh  (a,  a)  -  a,  ah  (a,  a) 
n'  1  '  1  n'  1  ’ 


(40a) 


This  can  be  compared  with  Rayleigh' s  result  for  a  fixed  rigid  sphere,  which  is 
obtained  from  the  boundary  condition 


(41) 


For  a  movable  sphere,  Equation  (40a)  is  still  valid  except  for  n  =  1.  The 
equation  for  A^  should  be  deduced  from  (33)  by  applying  the  same  limiting 
process  with  the  result 

(1  -»?)),(“,&)  -  «.aj  (o  a) 

A  =  -  3id>  - - — ± - ± £ ± _  / 

1  0  (1  -  q)h1(aia)  -  aiah2(aia)  ' 

Equation  (42)  can  be  converted  to  the  results  given  by  Lamb8  who  first  took  into 
consideration  the  rigid  translation  of  the  sphere. 
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For  viscous  fluid,  there  is  a  simple  analogy  between  waves  in  an  elastic 

37 

solid  and  a  viscous  fluid.  A  complete  discussion  has  been  given  by  Epstein. 

IX.  SCATTERING  CROSS  SECTION 

The  scattering  cross  section  is  defined  as  the  ratio  of  the  total  energy 

scattered  per  unit  time  by  the  obstacle  to  the  energy  per  unit  area  carried  per 

unit  time  by  the  incident  wave.  In  the  present  case,  the  total  energy  propagated 

1  2 

across  a  spherical  surface  of  radius  b,  concentric  with  the  scatterer  is  ’ 


u  *  +  T 
r 


rfl 


r*  u 
rr  r 


r*  u 
r6  0 


b  sin  0  ddd<P 


(43) 


The  energy  carried  by  the  incident  wave  per  unit  area  is 


q«>. 


iw 

2A 


\  \  (t  u*  -  r*  u  )  dxdy 
JJ  '  zz  2.  zz  z'  J 
A 


(44) 


The  superscript  (r)  and  (i)  of  the  stresses  and  displacements  in  (43),  (44) 
respectively  are  omitted,  and  an  asterisk  indicates  complex  conjugates.  With 
the  values  of  r  ,  u.  in  (13)  and  (14)  substituted  into  these  equations  the  scattering 
cross  section  is  found  to  be 


y  =  s!2 

(i) 

q 


=  4  7T 


2n  + 1 


n=0 


2  1  4> 

La!  0 


A  2 

IT  I  +  n(n+l) 


Ot  B 

r  i 


n 


(45) 


It  follows  then  y  has  a  dimension  of  an  area,  hence  the  name  "cross  section.  " 

By  comparing  the  coefficients  A^,  with  various  n' s  given  by  (22),  (28), 

(35)  and  (39),  it  is  seen  that  only  the  ones  of  n  =  0,  1,  2  have  the  same  order  of 

magnitude  for  small  Q^a.  with  different  kinds  of  spherical  obstacles  in  an 

elastic  solid,  the  scattering  cross  sections  all  can  be  expressed,  with  the  terms 

4 

of  order  higher  than  neglected,  as 
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y  =  (W9)g  a2 (a  a)4,  (46) 

In  which  g  takes  the  following  values: 

(1)  Elastic  Sphere  (\*  ^  p ^2*  P2  811  flnite) 


(2)  Fluid  Sphere  (p  =  n^/fx^  —  0) 
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(4)  Cavity  (p  —  0,  —  0) 


Since  q  is  a  dimensionless  quantity,  it  can  be  concluded  that  the  scattering 
cross  sections  are  all  proportional  to  the  sixth  power  of  the  radius  of  the  obstacle, 
and  inversely  proportional  to  the  fourth  power  of  the  incident  wave  length. 
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Figura  1 
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